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Abstract 

Using molecular dynamics simulations we study the be- 
havior of a dilute solution of strongly charged polyelec- 
trolytes in poor solvents, where we take counterions ex- 
plicitly into account. We focus on the chain conforma- 
tional properties under conditions where chain-chain in- 
teractions can be neglected, but the counterion concen- 
tration remains finite. We investigate the conformations 
with regard to the parameters chain length, Coulomb 
interaction strength, and solvent quality, and explore in 
which regime the competition between short range hy- 
drophobic interactions and long range Coulomb interac- 
tions leads to pearl-necklace like structures. We observe 
that large number and size fluctuations in the pearls and 
strings lead to only small direct signatures in experimen- 
tal observables like the single chain form factor. Further- 
more we do not observe the predicted first order collapse 
of the necklace into a globular structure when counterion 
condensation sets in. We will also show that the pearl- 
necklace regime is rather small for strongly charged poly- 
electrolytes at finite densities. Even small changes in the 
charge fraction of the chain can have a large impact on 
the conformation due to the delicate interplay between 
counterion distribution and chain conformation. 

1 Introduction 

Polyelectrolytes (PEs) are polymers which have the abil- 
ity to dissociate charges in polar solvents resulting in 
charged polymer chains (macroions) and small mobile 
counterions Because of their great relevance in tech- 
nical applications as well as in molecular biology they 
enjoy an increasing attention in the scientific commu- 
nity [21 El El El- The combination of macromolecular 
properties and long-range electrostatic interactions re- 
sults in an impressive variety of phenomena which makes 
these systems interesting from a fundamental point of 
view. 

In this paper we focus on the special case of poly- 
electrolytes under poor solvent conditions. The rea- 
son for this is that a large number of polymers are 
based on a hydrocarbon backbone for which water is 
a very poor solvent. The solubility in water is often 
only given due to their charged side groups. Impor- 
tant examples arc sulfonated poly-styrene (PSS), poly- 



methacrylic acid (PMA), DNA and virtually all proteins. 
The poor solvent conditions give rise to a competition 
between the attractive interactions due to the poor sol- 
ubility of the backbone and the electrostatic repulsion 
of the PE charges. This can lead to elongated strings 
of locally collapsed structures (pearls), commonly called 
pearl-necklaces. Such necklace conformations have been 
proposed on the basis of fluorescence studies [Hj. They 
have also been predicted in terms of scaling arguments in 
refs. 0121 EH for a weakly charged single chain PE. Later 
this has been extended to strongly charged chains and 
finite concentrations ^3 ^2 H'2I E3 ■ Also other theoret- 
ical approaches support the existence of necklace confor- 
mations ^] El EH] ■ The scaling approach of ref . [S] was 
supplemented with a Monte Carlo simulation of a sin- 
gle chain that shows a cascade from one to two to three 
globules with increasing strength of the electrostatic re- 
pulsion. In their study the polyelectrolyte was weakly 
charged and every monomer carried a fractional charge. 
The chain was studied in the infinite dilution limit, where 
counterions have not to be taken into account. The for- 
mation of the necklace structure is due to the Rayleigh 
instability of a charged droplet which leads to a split once 
a critical charge is reached. The size of the pearls is de- 
termined by the balance between electrostatic repulsion 
and surface tension. The distance between two pearls is 
governed by the balance of the electrostatic pearl-pearl 
repulsion and the surface tension. 

Some aspects of the theoretical pearl-necklace picture 
have been confirmed by simulations using the Debye- 
Huckel approximation |17II18| and with explicit counte- 
rions 

123122 E2 ESI El- However, there is up to now 
no clear experimental proof for the existence of necklace 
chains. Conformational chain properties have been ob- 
served which seem to be consistent with the necklace 

picture ElESESIEZlEiESlEni- 

Also the scaling of the peak position of the structure 
factor q* with the polymer density p has been thor- 
oughly investigated which however reveals only prop- 
erties that occur in the semi-dilute regime of interact- 
ing chains. For good solvent chains in the semi dilute 
regime the exponent [3 in the scaling relation q* oc p@ 
was measured to be (3 = 0.5 which is also the theo- 
retically predicted value whereas for poor solvent 
polyelectrolytes no single value for the exponent seems 
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to exist. The experimental values vary between j3 = 0.3 
and /3 = 0.5 (321 1331 15*^1 I35| and show a dependence on 
the charge fraction of monomers which is implicitly also 
responsible for the poor solvent parameter of the chain. 
The theoretical predictions show (for fixed poor solvent 
and fixed charge fraction) also a complicated transition 
from a /3 = g regime via a (3 = ^ into a crossover scal- 
ing |12l ll.'ij , whereas recent simulations measured a con- 
stant exponent j3 = g for the whole concentration range 
from semi-dilute up into the dense regime |24| , so its fair 
to say that things are far from being well understood. 
However, even for the dilute concentration regime, the 
situation can be more complicated as envisioned in the 
scaling approach, since the entropy of the chain and of 
the counterions as well as the electrostatic interaction 
between counterions and the PE charges have to be taken 
into account. In our previous shorter communications we 
reported large conformational fluctuations |22l 1231 l2"4"j in 
the pearl-necklace structures and showed that they were 
responsible for the absence of strong signatures in the 
force extension relation and in the form factor. 
The aim of the present simulational study is a more de- 
tailed investigation of the structure of strongly charged 
polyelectrolytes in poor solvent in the dilute concentra- 
tion regime where the chain-chain interaction is weak 
so that one deals effectively with single chain proper- 
ties. Our focus will be a thorough data analysis of the 
observed pearl-nccklacc conformations. To this end we 
had to develop a new cluster recognition algorithm that 
is capable to characterize these interesting conformations 
automatically from our simulated configurational data. 
Moreover we look at the stability of the pearl-necklace 
conformations in the presence of condensed counterions, 
and perform a study of the Coulomb induced collapse 
transition. We will compare our results with predictions 
from scaling theory and will discuss the validity range 
of the scaling approach for strongly charged chains at 
finite density. We then attempt to give a preliminary 
phase diagram for the systems studied. The last part 
will be devoted to some experimentally accessible ob- 
servables, like characteristic chain size ratios and form 
factor. This should be helpful in supporting the eval- 
uation of experimental data in terms of pearl-necklace 
signatures. 

Our paper is organized as follows: After explaining the 
used simulation method in sec. [21 and giving a short 
overview over the simulated systems in sec. [3] we will 
discuss our data analysis methods in sec. 0] In sec. |S] 
we will compare our results in the pearl-necklace regime 
with predictions made by scaling theories. In the next 
sec. we will quantify the role of fluctuations, and de- 
scribe in sec. [7|our view of the Coulomb induced collapse 
transition and a preliminary phase diagram for the range 
of our simulation parameters. We then discuss in detail 
some measurable observables in sec. |H1 and finally end 
with our conclusion in sec. El In the appendix we give a 
detailed overview of the simulated systems, including a 



list of parameters and results for selected basic observ- 
ables. 

2 Simulation method 

Our model of a PE solution and our molecular dynam- 
ics approach has been described previously in detail 
in ^3 123 It consists of Np flexible bead-spring 

chains with N m monomers and N c counterions which arc 
located in a cubic simulation box of length L with peri- 
odic boundary conditions. A fraction / of the monomers 
is monovalcntly charged (y m = 1). Thus the total charge 
per chain is Qp = fN m . The number of counteri- 
ons which are also monovalently charged [y c = —1) is 
chosen such that the overall system is electrically neu- 
tral. Densities are given either as monomer density 
Pm = N P N m /L 3 or charge density p c = 2 fN P N m /L 3 . 
All particles interact via a Lennard-Jones (LJ) potential 
^ e i(f) 12 ~ (f) 6 — c ] f° r distances r < R c and zero else- 
where. The constant c is chosen such that the potential 
value is zero at the cutoff R c , and e is a measure of the 
solvent quality. Monomers interact up to R c = 2.5c giv- 
ing them a short range attraction which can be tuned 
by changing the value of e. The counterions interact via 
a purely repulsive LJ interaction with R c = 2 1 / 6 er. The 
units of length, energy and time are a, e, and r, respec- 
tively. 

The chain monomers are in addition connected along the 
chain by the finite extendible nonlinear elastic (FENE) 
bond potential of the form -14fc B Tln[l - (r/(2er)) 2 ] 
which results in an average bond length few 1.1 a. 
Charged particles with charges qi and qj at separation 
Tij interact via the Coulomb energy k^Ti^qiqj jr^ where 
the Bjerrum length is defined as £b = e 2 / (A-Kesiok-QT) 
(e: unit charge, eo and eg: permittivity of the vacuum 
and of the solvent). The Coulomb interaction was cal- 
culated with the P3M-algorithm [SZIEH1, tuned to force 
accuracies which are much higher than the thermal noise 
level. 

A velocity Verlet algorithm with a standard Langevin 
thermostat is used to integrate the equation of mo- 
tion |2H1 (friction coefficient T = t -1 , time step At = 
0.0125t). Thus the solvent is only implicitly present via 
its permittivity £5, the friction constant T and the sol- 
vent quality parameter e in the LJ potential. 
The simulation time after equilibration for all systems 
was at least 100 times the measured correlation time for 
the end-to-end distance Re and the chains centers of 
mass diffused at least several radii of gyration Rq. The 
osmotic pressure p was measured to be always positive 
and additional simulations over a large density range |2"1] 
showed that the pV diagram is convex at all densities, 
thus our simulations are stable, reach true thermal equi- 
librium, and reside in a one phase region. The volume 
density inside the pearls does not exceed 0.47 which is 
below the glass transition. We therefore are certain that 
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also the pearl formation and restructuring was observed 
in equilibrium. 

3 Simulated systems 

We have simulated our poor solvent polyelectrolytes 
mainly as a function of the parameters chain length N m , 
solvent quality e, strength of the electrostatic interac- 
tion £b and the charge fraction /. All simulations are 
performed in dilute solution such that the interaction 
between the chains is small. For the system with the 
longest chains (N m — 478) we have a chain extension 
i?E ~ 60cr and a chain-chain separation r cc w 252a which 
was calculated for a random packing of spheres according 
to 
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The screened renormalized monopole interaction C/^ H 
between two chains can be estimated for this case to 
be of the order fceT. For our estimate we have used a 
crude approximation on the Debyc-Hiickel level. The ef- 
fective chain charge Qp, c ff ~ 64 is calculated by using 
the counterions within a shell of 3<r around the polyelec- 
trolyte to renormalize the bare charge Qp = 160. The 
density of the free counterions p c leads a screening con- 
stant k = \/Att£bPc ~ 0.009er~ 1 and as the interaction 
potential we take C/° H = 4s<3p,eff e ~ Krcc r cc ~ 2M B T. 
For this estimate one has to keep in mind that both the 
used practical definition of the effective charge as well 
as the screening concept can not be founded on physical 
principles. 

Details about the used parameters and measurements of 
some basic observables for all simulated systems can be 
found in appendix A. We have grouped the simulations 
into series depending on the investigated parameters. 
The chain length dependence is studied in simulation 
series which differ in the line charge density /. For series 
Al we use a charge fraction / = i , and the chain length 
is varied over one decade in steps of 48 from N m = 48 to 
478. Scries A2 and A3 are performed at N m = 100, 200 
and 300 with / = i at slightly different densities. 
To study the dependence on the solvent quality we have 
varied the short range attraction by changing the LJ pa- 
rameter e between 0.0/cbT and 2.0fce2 n for medium sized 
chains of length N m = 238 (series B). The O-point for 
this model was determined to be at e(Q) = 0.34fceT |19j . 
For most of the other simulations (series Al, A2, A3, CI, 
C2 and C3) we used e = 1.75fceT which is thus deep in 
the poor solvent regime. This value is chosen for practi- 
cal reason. In this regime we have found relatively large 
and stable pearl-necklace conformations which are more 
easy to investigate. 

The effect of the Coulomb interaction is studied via 
changing the Bjerrum length £b and the charge fraction 
/. We simulated chains with length N m = 199 in three 



series CI, C2 and C3 with charge fractions / = 1, / = | 
and / = 3 respectively. They are performed at the same 
charge density p c = 5x 10~ 5 er~ 3 but they contain differ- 
ent numbers of counterions N c = fN m Np corresponding 
to the number of charges on the PEs. All three series 
start out at £b = Oct. In Series CI Ib ranges up to 10cr, 
in series C2 and C3 up to 9a. In addition there is a 
simulation set with shorter chains N m = 94 with f = k 
and a different solvent parameter e = LSfceT 1 which we 
only use for the phase diagram (series C4). 
The last series contains four simulations with the same 
value for the scaling variable £b& -1 / 2 = 0.25 but differ- 
ent / = 1, J! I an d \ (series D). 

4 Observables and data analysis 

In this section we define our measured observables and 
explain in detail how we analyzed our simulated PE con- 
formations. Especially we present the cluster recognition 
algorithm which we used to automatically classify differ- 
ent pearl-necklace structures. 

We denote the position of monomer i with r^ and the 



distance between two particles i and j with 



The 



center of mass for the chain is then R s = jj— X^=!i r * 
and the center of mass coordinates are = r.; — R s . For 
the chain extension we use the end-to-end distance 



i?| = (n - r Nl[L ) , 



the radius of gyration 
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and the inverse hydrodynamic radius 



1 Nm 1 

D-l = 1 V — 

N m (N m - 1) ^rf n 



(2) 



(3) 



(4) 



i j 



Note that this definition corresponds to the short time 
diffusion behavior of polymers. For more information on 
this topic see ref. |4()j . 

For a first structure classification we use two charac- 
teristic ratios between the different chain extension ob- 
servables. The first characteristic ratio is defined as 
r = (Re/Rg) 2 - The second characteristic ratio is 
a = Rq/Rh which has the advantage of being exper- 
imentally accessible |41j . 

We will also compute the spherically averaged form fac- 
tor, sometimes also called single chain structure factor, 
Si(g) which can be measured in scattering experiments: 



N m N m 



i=l j = l 



sin(qrj 



(5) 



To describe the counterion distribution around poly- 
electrolytes we use an integrated counterion distribution 
P(r). The distance di of a counterion i from a chain is 
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defined as the distance of the counterion from its closest 
monomer in space. We denote the set of counterions be- 
longing to one chain with C. From this we can calculate 
P(r) as 

P(r) = J- f dr'^Sir'-di) (6) 

P(r) denotes the fraction of counterions which are in- 
side a shell with radius r around a polyelectrolyte chain. 
This definition can be used for a large variety of chain 
conformations. 

For an automated analysis of all types of pearl-necklace 
structures appearing in our configurational data we need 
a tool to determine for each monomer, to which pearl or 
string it belongs, and what is the total number of such 
substructures. First we should state that there is no 
sharp definition of a pearl or a string. So what we need is 
a practical approach to the problem. Our guideline in the 
development of an automated tool for the identification 
of pearls and strings is that the result should be close 
to the result which would be obtained by looking at the 
conformation by eye. 

Before we give a detailed explanation of the used clus- 
ter recognition algorithm we want to present some other 
methods that are based on observables which could be 
accessible with experimental methods. In this way one 
can also judge how easily theses structures can be ob- 
served by current experimental techniques. 
The local monomer concentration p m (r) of a chain is 
given by p m (r) = J2i m ^( r ~ r i)- To distinguish the 
different substructures, pearls and strings, it is better to 
use a coarse grained local monomer concentration Pm( r ) 
which is defined as 

/4(r) = ^ / drV m (r'), (7) 

where V c is a spherical volume with radius r c around 
r. This observable is larger in a pearl than in a string. 
One can calculate Pm( r ) along the backbone as a func- 
tion of the monomer positions Yj. This is done for a 
chain with length N m = 382 from series Al whose snap- 
shot is shown in fig. The resulting coarse grained 
local monomer concentration with r c = 4a is shown for 
comparison directly underneath in fig. 
A second approach uses local distances between the 
monomers which could be probed by NMR experi- 
ments |30| . The local distance between monomer j and 
the monomer j + n which is n monomers apart on the 
backbone is defined as 

r n(j) = \J (rj - r j+n ) 2 . (8) 

From scaling arguments one expects that the local dis- 
tances inside the compact pearls scale with r„ ~ n 1 / 2 
for small n whereas in the extended strings it scales with 
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Figure 1. a) Snapshot of a polyelectrolyte chain 
with iV m = 382 from series Al. b) Coarse grained local 
monomer concentration p^C?) with r c = 4a along the 
backbone, c) Local distances r n (j) for n = 10 along 
the backbone, d) Analysis with the cluster recognition 
algorithm. Monomers in pearls are grey and monomers 
in strings black. 



r n ~ n. The symbol ~ is used to state a scaling rela- 
tion which ignores numerical prefactors. For the con- 
formation shown in fig. ^a, the local distances along the 
backbone are shown for n = 10 in fig. 
In both cases one can see the position of the 4 pearls 
and 3 strings. The problem is that the variation of both 
p^(r) and r n (j) is of the same order as the difference 
between the mean values in pearls and strings. 
Since the two previous methods have several problems 
with the structure recognition of pearl-necklace confor- 
mations we have developed a more reliable algorithm 
based on a simple cluster recognition. The question 
whether a group of monomers forms a cluster is often 
connected to a distance criterion, e.g. monomers with a 
distance smaller than a critical value belong to the same 
cluster. For a polymer one has also to take into account 
the chain connectivity. This implies that a pure distance 
criterion is not sufficient. Therefore we require in addi- 
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tion that there are a certain number of bonds between a 
pair of monomers along the chain contour. The resulting 
algorithm is iterative and contains the following steps: 

1 . At the beginning every monomer is a cluster of size 1 
(size = number of monomers belonging to a cluster). 

2. Two clusters CI and C2 are merged if they contain 
a pair of monomers ij with i £ CI, j G C2 and 
fij < T c - In addition i and j are further apart than 
n c bonds along the chain contour: \i — j\ > n c . 

3. Step 2 is repeated between all clusters as long as 
one finds clusters that have to be merged. 

4. Remove loops: a cluster CI where all monomers are 
inside (along the chain contour) of another cluster 
C2 are merged. Note that this step is only suitable 
for polyelectrolytes, but e.g. is not applicable for 
p olyampholy t es . 

5. Practical definition of pearls: all clusters with a size 
larger than or equal p c are pearls. Pearls which 
are connected directly along the chain contour are 
merged. 

6. Practical definition of strings: all clusters with a 
size smaller than p c belong to strings. Strings which 
are connected directly along the chain contour are 
merged. 

7. Remove dangling ends and merge them to the end 
pearls. We do this because we have so far not seen 
dangling ends containing more than 3 monomers 
which could be seen as an extra string. 

Thus the algorithm has in principle three free parame- 
ters: r c , n c and p c . But looking at the involved struc- 
tures, one can establish relations between these parame- 
ters which can be used as a rough guide in choosing them. 
From scaling arguments wc know that the distance of two 
monomers i and j with a distance \i — j\ along the chain 
contour scales as \i— j\ x / 2 inside the pearls and with 

1/2 

inside the strings. Thus we can choose bn c < r c < bn c 
to distinguish the two cases. For a weakly charged chain 
one can determine a suitable value for p c with help of the 
pearl size defined in eq. 1131 But this is not possible for 
our strongly charged systems with their subtle depen- 
dence on the counterion distribution. The whole data 
analysis in this paper is done with an empirical parame- 
ter set: r c = 2.1a, n c = 6 and p c = 9. Beside extended 
visual checks we have tested that small changes of the 
three parameters do not have a significant effect on the 
final result. The derived sizes for the substructures con- 
tain a systematic error of ±4 monomers coming from 
the two outer monomers on each side of a substructure, 
for which one can not decide whether they belong to 
the next pearl or the next string. A typical result for a 
structure type with 4 pearls is shown in fig^d. The four 



pearls and three strings contain the following numbers of 
monomers: 90 -8- 94 -6- 77 -9 - 98. We do not claim 
that this is the fastest or the best way to identify pearl- 
necklace structures, but it worked well when compared 
with visual checks and thus served our purpose. For an 
average pearl size larger than 30 monomers it yields a 
reliability well above 95%. 

5 Scaling 

One of our goals in this paper is to show to what extend 
scaling theories that are made for long chains at infinite 
dilution can be expected to work for dilute PE systems 
with finite length at finite density. The scaling theory 
predicts the dependencies of observables like the chain 
extension on various parameters, e.g. N m , £b, f and 
the reduced temperature r r . Here we give only a short 
overview of some results of the scaling theory for the 
pearl-necklace regime of polyelectrolytes. For the pearl- 
necklace regime one finds the following relations EH 
HT] : End-to-end distance Re 

Rv^N^/H^fr-^, (9) 
number of pearls np 

np-N^-Hnfr- 1 , (10) 

string length Is 

h^lH^I 2 r X r}l\ (11) 
pearl size (radius) rp 

r P ~ 6 4 / V /3 r 2/3 , (12) 

pearl size (number of monomers) gp 

g P ~ b£^f- 2 T T , (13) 

density inside the pearls pp: 

P p ~ 6- 3 r r . (14) 

Note that rp is independent of the solvent quality, and 
that our solvent quality parameter e is proportional to 
the second virial coefficient of the LJ potential and hence 
also proportional to the reduced temperature r r . 

5.1 Scaling variable chain length N m 

The linear scaling of Re with N m is caused by the elec- 
trostatic repulsion of the chain charge and is dominated 
by the string length Zg. Due to the finite length of our 
systems we have to correct the scaling relation given in 
cq.[3] For small chain length, JV m < .gp, the chain confor- 
mation consists of one pearl. The size of this pearl scales 

1 /3 

as Nm ■ Thus we have to replace N m by (N m — gp) in 
cq.EO In fig. [51a linear relation between Re and (N m —gp) 
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100 200 300 400 

N m"9p 

Figure 2. Scaling plot for i?E versus (N m — gp) (see 
eq. |3J. The series Al, A2 and A3 differ in the charge 
fraction / and the density p c . One can see that Re 
scales linearly with (N m — gp), but the prefactor is not 
constant. 



can be observed. But the full scaling relation from eq. |U] 
including the parameters specifying the Coulomb inter- 
action, namely /, is not valid. Series Al is performed at 
/ = 3 and its density is in between that of series A2 and 
A3 which have both / = |. 

We find a similar result for the number of pearls np in 
the pearl-necklace regime. The scaling relation for np 
versus (N m — gp) (see ea. ll()(l is shown for the data from 
series Al, A2 and A3 in fig. Again there is a linear 
relation between np and N m but different prefactors for 
different values of / and p c . 

To get some visual impression, we show in fig. 0] some 
snapshots of simulations with different chain length N m 
from series Al. The discreteness of the pearl number np 
does not play a significant role for sufficiently long chains 
which give rise at least to a dumbbell. Then our data 
is in accord with the scaling relations for the quantities 
pearl size, string length and pearl density, namely they 
are constant within the statistical error. We find for 
the systems of series Al with N m > 200: gp = 78 ± 4, 
p P = 0.67 ± 0.04cr~ 3 and l s = 7.3 ± 2cr. As wc will show 
later, the discreteness of the number of pearls is smeared 
out due to fluctuations between different structure types 
in a way that the average quantities can maintain their 
optimal values, compare also sectional 

5.2 Solvent quality e 

In the simulation series B we have tested the behavior of 
polyelectrolyte chains upon changing the solvent quality 
via the Lcnnard-Jones parameter e. Since the practically 
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Figure 3. Scaling plot for np versus (N m — gp) (see 
cq.rTul for series Al, A2 and A3. 




N m =286 . 

e * • 




Figure 4. Snapshots from series Al with different 
chain length N m and different number of pearls np. 



usable range for this parameter is small, ranging from 
e = OfcpT to 2.0fcBT, it is not possible to test the scaling 
predictions over a large parameter range. In addition the 
range of e values for which we observed pearl-necklace 
structures is even smaller, namely from e = l.OfcpT 1 to 
2.0fcp,r. For values of e < l.OksT one leaves the poor 
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solvent regime and for values e > 2kp,T we encounter 
simulational problems, e.g. kinctically frozen states. 
The dependence of Re and np on e is shown in fig. [^1 and 
n 8- El F° r both observables the scaling predictions seem 
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Figure 5. Dependence of Re on the solvent quality 
parameter e (simulation series B). 
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Figure 6. Dependence of np on the solvent quality 
parameter e (simulation series B). 



to hold for a small e-region. Re scales linear with e -1 / 2 
in the range from e = 1.25kp,T to e = Llbk^T (see eq.EJ. 
The regime where np scales as e" 1 is even smaller and 
extends from e = \A5k B T to e = 1.75k B T (sec ea.Hn|). 
In the same regime also the density inside the pearls pp 
which is plotted against e in fig. shows a linear de- 
pendence on e as predicted by cq. 1141 However, a more 
reliable conclusion may be drawn from the scaling rela- 
tion for the pearl size rp, since rp should be independent 
of the solvent quality (see ea. H2f) . We therefore have cal- 

o 1/3 

culated the pearl size as rp = jirgppp from our sim- 
ulation data and plotted it versus e. As can be seen in 
fig'-El^p increases monotonically with e, hence shows an 
unexpected dependency on e. We are lead therefore to 




1.5 1.6 
£[(k B "0] 

Figure 7. Dependence of pp on the solvent quality 
parameter e (simulation scries B). 
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Figure 8. Dependence of rp on the solvent quality 
parameter e (simulation scries B). 



conclude that the scaling predictions for the dependency 
on the solvent quality do not work here. They fail for 
the most definite test case, the independence of the pearl 
radius rp of the solvent quality. Nevertheless in a small 
regime for e/(fceT) = 1.45... 1.75 they seem to work 
approximately for some observables due to a fortuitous 
error cancellation. 

To elucidate the role of the finite counterion density 
around the chains we calculated the integrated ion dis- 
tribution P(r) (see eq. EJ). In fig. P(r) is shown 
for different solvent qualities. The counterion distri- 
bution changes greatly close to the chain with e. In 
the pearl-necklace regime between e/(ksT) = 1.0 and 
e/(fceT) = 1.85 which is the same regime where we find 
also the strongest change in the chain extension we ob- 
serve that the counterion distribution is changing the 
most. The fraction of counterions being very close to the 
chain which hence can be called condensed is also vary- 
ing strongly. Only the top four curves show an inflection 
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Figure 9. Integrated counterion distribution P(r) 
around polyelectrolyte chains subject to different solvent 
qualities. 



point which is a sign of counterion condensation [421 143) 
(see also sec. 16.10 . Scaling theory assumes usually that 
the counterion condensation depends only on the Man- 
ning parameter £ = ^b^ -1 / an d is hence supposed to 
be independent of the actual conformation of the chain. 
Note that the total series B is performed at £ = 0.5 and 
is thus expected to show no counterion condensation at 
all according to the standard Manning-Oosawa ^] 14*5*] 
concept which requires £ > 1. 

5.3 The Coulomb parameters: £& and / 

As we demonstrated in the previous subsection 15.21 the 
interplay between counterions and chain conformation 
influences greatly the behavior of PEs in poor solvent. 
Therefore also a thorough investigation of the parame- 
ters determining the Coulomb interactions is necessary. 
This section treats the simulation series CI, C2 and C3 
where we have investigated the £b dependence but us- 
ing different charge fractions / = 1(C1), / = 1/2(C2) 
and / = 1/3(C3). All simulations are performed at the 
same charge density p c . We remark that the valences 
of the charged particles is yet another important inde- 
pendent parameter. But since this is also not included 
in the scaling picture and would even further compli- 
cate the picture we will leave this for another study. We 
will also look at the behavior of these systems with re- 
gard to parameters combined of £b and /, namely the 
Manning parameter £ and the scaling variable £b& -1 / 2 - 
The Manning parameter is important for the electro- 
static field that counterions would experience around a 
long stretched chain and its value determines in a first 
approximation the onset of counterion condensation. In 
the framework of scaling theory simulations performed 
at the same value of ^b^ -1 / 2 should be identical. Nev- 
ertheless our simulations show a big difference between 
the three series. The dependence of Re on £b is shown in 
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Figure 10. Dependence of Re on the Bjerrum length 
£b for the series CI, C2 and C3. The inlay shows the 
dependence of Re on the scaling variable ^b& -1 / 2 - 



ues of £b the counterions and their interaction with the 
chain conformation do not play a significant role. The 
prediction from scaling theory (see eq.|HJ) is that our data 
should collapse on a single master curve. However, our 
data show that this is only true if the scaling variable 
£b^ 1 f 2 is smaller than 0.1. This can be seen in the in- 
set in the upper right corner of fig. 1101 where we plotted 
the same data versus £bI>~ 1 f 2 . Upon a further increase 
of £Bb 1 f 2 the values of the end-to-end distances diverge 
rapidly. Whereas the chains of series C3 (/ = 1/3) al- 
ready start to shrink at ^b& -1 / 2 = 0.15 the chains of the 
other series still expand. 
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Figure 11. Dependence of np on the Bjerrum length 
£b for the series CI, C2 and C3. 



The different values of / can be seen as different schemes 
for the discretization of the backbone charge. This 
has an effect on the correlations between the charges 
which influences the chain conformation and thus also 
the counterion distribution. In fig. 1121 we have therefore 
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plotted the integrated counterion distribution P(r) for 
several values of the Manning parameter £ for series C2 
and C3. Already the curves for £ = 0.5 which corre- 
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Figure 12. Integrated counterion distribution 

P(r) around polyelectrolyte chains for different ^-values. 
Open symbols show data from series C2 with / = 1/2 
and filled symbols from C3 with / = 1/3. 



spond to £b = If for series C2 and £b — l-5<7 for series 
C3 show a pronounced difference in P(r). Even though 
the effect of / on the charge charge correlations may be 
small it is enhanced strongly by the interplay between 
the chain conformation and the counterion distribution. 
When the counterions only slightly move towards the 
chain due to correlation effects, the effective charge of 
the chain will shrink and so does the end-to-end distance. 
A higher effective Manning parameter £r := Qp&/i?E 
follows which again attracts more counterions towards 
the chain. The same holds for the opposite way. This is 
the same mechanism which we have already seen for the 
dependence on e in the previous section. The difference 
in P(r) becomes smaller with increasing £ which is due 
to the gradual collapse of the polyelectrolyte chain, since 
then in all cases most of the counterions are close or even 
inside the chains. 

In series D we have performed four simulations with 
^B& -1 / 2 = 0.25, but different values for £b and /. In 
fig. El w e show the behavior of np as a function of £. 
Decreasing / and thus increasing £b has again a drastic 
effect on the chain conformation leading to collapse of 
the chains for large £b and small /. Even though the 
number of pearls for the two first points with / = 1 
and / = 1/2 do not differ much the chain extension 
shows a large difference, namely Re — 59<r for / = 1 
and Re = 32er for / = 1/2. This is also reflected in a 
large difference of the pearl sizes and the string lengths. 
The difference between the systems can again be traced 
to a quite different distribution of the counterions as can 
be inspected in fig. 1141 where we show P(r) for three of 
the simulations of series D. It is really surprising that 
such a small change, basically in the discretization of 




Figure 13. Dependence of np on the Manning para- 
meter £ for the simulation scries D with constant scaling 
parameter £b& -1 / 2 — 0.25. 
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Figure 14. Integrated counterion distributions P(r) 
are shown for different charge fractions / with constant 
£ B b- 1 f 2 = 0.25 (series D). 



the charges, has such a big impact on the chain con- 
formations. This again shows that the delicate balance 
between repulsive and attractive forces is very sensitive 
to subtle changes. These effects are definitely not cap- 
tured by the parameter ^b^" 1 / 2 which is used in scaling 
theories. We will come back to the collapse discussion in 
scc.0where we will also give an overview of the occurring 
conformations. 



6 Fluctuations 

The data analysis shows that pearl-necklace conforma- 
tions are very soft objects which display large fluctua- 
tions on all length scales. We will measure in detail the 
extent of the various fluctuations, since this is important 
for the interpretation of experimental measurements on 
polyelectrolyte solutions (see section [HJ. 
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6.1 Fluctuations of the structure type - Coexis- 
tence 

In most of our simulated systems we find coexistence of 
pear-necklaces with different number of pearls which we 
call structure types. It is important to note that the ob- 
served coexistence is not caused by frozen meta-stable 
states. We have excluded this possibility by observing 
individual chains over a larger period of time each show- 
ing a large number of transitions between different struc- 
ture types during the simulation time. The typical time 
evolution of the structure type np and the radius of gy- 
ration Rq for an individual chain is cxemplarily shown 
in fig. 1151 In this case the chain mainly fluctuates be- 
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Figure 15. Time evolution of the structure type (num- 
ber of pearls) and the radius of gyration Rq of a single 
chain (N m ~ 382, series Al). 



tween structures with 4 and 5 pearls, but we also see 
a significant fraction of structures with 3 and 6 pearls. 
Looking at Rq one can see a certain correlation with the 
structure type which is reflected also in the mean value 
for different structure types: Rq* = (14.7±1.3)er, Rq* = 

(16.1 ± 1.5)cr, R ( q } = (17.5 ± 1.6)cr, R ( £ } = (19.5 ± 1.7) a. 
Here and in the following the superscript denotes that 
a observable is measured only for conformations with a 
certain number of pearls np. But nevertheless such a 
simple chain observable is not suited for structural dis- 
crimination as one can see from the given mean devia- 
tions. The same holds for other chain obscrvables like 
i?H, the characteristic ratios r and a. 
One could be tempted to explain the coexistence of two 
structure types by a simple finite size argument, namely 
that the ratio of N m and the optimal pearl size gp is 
not always an integer. In fig. 1161 we plot the probability 
distributions P(np) of the structure types for three dif- 
ferent chain lengths. We clearly observe a coexistence 
range containing up to four different structure types. 
This means that the different observed structure types 
can only have small differences in the free energy of the 
order of the thermal fluctuation spectrum. To confirm 



Figure 16. Probability distribution for the structure 
types found in systems with different chain length from 
scries Al. 



this we calculate from the shown probability distribution 
the free energy differences between structure types with 
n and m pearls according to the Boltzmann factor: 
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For the case with N m = 430 we find: AJ^ 45 ) = 
-1.32 feT, A^ 56 ) = 0.66 k B T, A^ 67 ) = 1.9k B T. All 
values are of the order of k B T which is consistent with 
our observed large coexistence regime. 
As we will argue these small differences in the free en- 
ergy between different structure types are mainly due to 
the interplay of the counterion distribution and the chain 
conformation. To elucidate the role of the countcrions 
we look at the counterion distribution around the PEs 
analyzed for each structure type separately. In fig. 1171 
we have plotted the integrated counterion distributions 
P(r) for different structure types for the system with 
N m = 430 from series Al. One observes more counteri- 
ons to be close to the chains with a smaller pearl num- 
ber. This can be easily understood if one looks at the 
far electrostatic field of the chain. For distances larger 
than the pearl-pearl separation rpp the chain can be seen 
as a charged cylinder with an effective Manning charge 
parameter £r„ := Qpb/Rp,. The end-to-end distance 



increases with increasing pearl number: R 



(4) 



48.8er, 



R^ = 54.1(7, R y °> = 60.1a and R^> = 64.0cr, hence the 
shorter chains have a larger effective Manning charge 



? (f>> 



?(7) 



parameter, namely £ 



(4) 



2-95, £g 



and d 7) 



2.67, $1 = 2.40 



2.25. 



Since the PEs in this regime are elongated structures 
and they carry a large effective line charge it might be 
worthwhile to compare the ion distribution with pre- 
dictions from Poisson-Boltzmann (PB) theory. This is 
also supported by the functional form of the integrated 
counterion distributions in fig. I17l which looks very sim- 
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Figure 17. Integrated counterfoil distributions P(r) 
for different structure types for the system with N m = 
430 from series Al. The arrow pointing to (17a, 0.53) 
shows the position of the inflection point for the struc- 
ture type with np = 4. 



ilar to that of an infinitely charged rod within the cell 
model jJB] . In the framework of the PB theory applied 
to the cell model for an infinitely long charged cylin- 
der 07j one can calculate the fraction f c of Manning- 
Oosawa condensed counterions as / c = 1 — |. 

In the simulation the fraction of condensed counterions 
can be read of as the value of P(r) where the P(r) curve, 
plotted as a function of log(r), has an inflection point. 
Using £ = £r e we find for the np — 4 structure that 
the cell model prediction for f c = 0.66 is higher than 
the value which can be read off the inflection point crite- 
rion, that is f c = 0.53. A perfect quantitative agreement 
can however not be expected since our system is finite, 
and has as such a smaller electric field, and secondly the 
pearl-necklace structures do alter the near electrostatic 
field of the chain and lead to an inhomogeneous counter- 
ion distribution. As a word of caution we remark that 
the bare value for the Manning parameter of the chain in 
Fig. ll7l is 0.5 which would not lead to any counterion con- 
densation at all in the simple Manning picture, and the 
functional form of an infinite rod with that charge value 
would simply not show any inflection point. This again 
demonstrates that poor solvent chains, due to the short 
range attraction between monomers, feature a higher ef- 
fective charge density. Since the Poisson-Boltzmann cell 
model can at least qualitatively explain the change in the 
counterion distributions a more refined version of this 
model seems to be necessary. A possible starting point 
could be to combine the PB rod and sphere geometry as 
has been suggested in ref. |48j . 

When more counterions are near the chain we find in 
terms of a charge renormalization that the effective line 
charge density of the chain decreases. This allows the 
chain to contract further which itself induces a stronger 
counterion attraction. Again a lower effective line charge 



density also increases the optimal pearl size. This re- 
sults in shorter chains with fewer pearls. The opposite, 
longer chains with more pearls, holds when the counteri- 
ons move away from the chain. This lowers the difference 
in the free energy between the different structure types, 
as was suggested above. 

Scaling theories have predicted a collapse of the pearl- 
necklace structure into a globule as soon as counterion 
condensation starts [491 1101 l5TJ| due to an avalanche 
behavior of condensing counterions that contracts the 
chain. In the investigated parameter range where we 
observe stable pearl-necklace conformations in our simu- 
lation data we did not see such a collapse transition. In 
the light of a recent study |51| this is not surprising, since 
it can be shown that this collapse depends on how easily 
the counterions can enter the globule, the ion concen- 
tration, and the strength of the electrostatic interaction. 
Also the necessary amount of condensed counterions to 
induce the collapse can be quite high. 
We also suggest that the fluctuations due to the pres- 
ence of the counterions lower the energy barrier between 
the different structure types. To confirm this sugges- 
tion one would need to analyze the transition frequency 
between the different structure types for chains with 
explicit counterions and a corresponding Debye-Hiickel 
chain with the same np. This is however outside the 
scope of the present investigations. 

6.2 Fluctuations of the substructures 

In this section we discuss fluctuations of the chain con- 
formation on a smaller length scale, namely the radius 
of the pearls rp and the pearl-pearl distance rpp. The 
probability distributions for the pearl size gp is shown 
in fig. 1181 and for the pearl-pearl distance rpp in fig. 1191 
Both distributions are quite broad. The relative stan- 
dard deviations for the distributions of all chains are 
Sgp = 0.32 and Sr PP = 0.22. 

As we have already seen for the overall chain extension 
also the size of the substructures arc influenced by the 
counterion distribution. The pearl size is decreasing with 
increasing number of pearls and increasing chain exten- 
sion as one can see from the arrows in fig. 1181 Still an 
explanation in terms of a charge renormalization is not 
sufficient for an understanding. This becomes more clear 
when we look at the number of counterions inside a shell 
around the chains. Rounded to integers we find for a 
shell radius of 3cr on average 56 counterions. Averaged 
separately for the different structure types we find 60, 
57, 53 and 50 counterions for structures with 4, 5, 6 and 
7 pearls. It is interesting to note, that the difference in 
this number of counterions is roughly constant for shells 
with radii ranging from 2a up to 30tr and is thus in- 
dependent of the exact definition of counterions called 
"condensed" to the chains. Thus we now use a practi- 
cal definition of an effective charge using the counterions 
which are closer than 3er to the next monomer. For the 
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Figure 18. Probability distributions P for the pearl 
size gp for the system with chain length N m = 382 from 
series Al. Shown is the distribution for all chains as 
well as the distributions for the different structure types. 
The arrows mark the mean value of the corresponding 
probability distribution. 
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Figure 19. Probability distributions P for the pearl- 
pearl distance rpp for the system with chain length 
N m = 382 from series Al. Shown is the distribution 
for all chains as well as the distributions for the different 
structure types. The arrows mark the mean value of the 
corresponding probability distribution. 



example this yields an effective pearl charge of 18, 15, 13 
and 11 charges for structures with 4, 5, 6 and 7 pearls 
respectively. Since the pearl size is only slowly varying 
one gets large differences in the electrostatic self energy 
of the pearls of different structure types. 
In the scaling Ansatz the distance between neighbor- 
ing pearls rpp is determined by the balance between 
the electrostatic repulsion between the chains and the 
energy one needs to pull monomers out of the pearls. 
The second term is related to the surface tension and is 
roughly constant. The term connected to electrostatics 



is more complex. We approximate it by the electrostatic 
energy of two point charges carrying the average pearl 
charge, separated by the average pearl-pearl distance. 
?*pp ranges from 15.4er for 4-pcarl structures to ll.lcr for 
7-pearl structures. Together with the above calculated 
effective pearl charges we can estimate the electrostatic 
repulsion between neighboring pearls to be 3l.9kBT for 
4-pearl structures and 16.7kBT for 7-pearl structures. 
Note that we neglect the energetic contribution of the 
other pearls. This difference can also not be explained 
by a standard Debye-Huckel potential since the Debye 
screening length, calculated from the bulk charge den- 
sity p c = 1 x 10 -5 er -3 in the system, is = 72.8<r and 
is thus much larger than the values for rpp. Of course, 
closer to the chain the counterion density is much larger 
than in the bulk which could be a hint to explain the 
observed differences. The counterion density is rapidly 
varying with the distance from the chain starting with 
w xl0 _1 <7 _3 at a distance of 1.5a and dropping below 
the bulk density at a distance of 50o\ 
The electrostatic self energy of the pearls and the elec- 
trostatic pearl-pearl repulsion show explicitly the dis- 
crepancies between a mean field approach and our sim- 
ulational results. It also clearly demonstrates that it 
is necessary to include the counterion distribution into 
the model. Moreover correlations between the counte- 
rions and the chain charges seem to play an important 
role. For instance, it can be shown that the counterf- 
oils preferentially accumulate between the pearls. Since 
the structure of polyelectrolytes in poor solvent turns 
out to be extremely sensitive to the inhomogeneities of 
the counterion distribution it should be used as a test 
case for the development of theoretical approaches be- 
yond the mean field level. A more detailed discussion of 
these inhomogeneities can be found in ref . |21l 15*2*] 

7 The sausage regime and phase space 

In order to get an overview of the parameter regimes our 
different simulation series are scanning through we have 
depicted them in a phase plot for polyelectrolytes using 
the parameters e and £ B ~ 1 as it was done for example 
in Ref. ^UJ- Some of the simulation series are marked 
with dashed lines in fig. [201 Close to the lines we have 
put some snapshots of the configurations to get a visual 
impression. 

Looking at series C3 crossing the entire ^B-range we 
find, starting on the right side at the neutral case with 
£b = O.Oer, a neutral globule. The chain is collapsed 
due to the poor solvent. With increasing £b the chain is 
getting charged, the self-energy grows, and we can ob- 
serve the Rayleigh instability at £b = 0.75ct where the 
globule splits into a dumbbell. At this point we find al- 
ready some counterions close to the chain, hence there 
are condensed counterions present. Further increase of 
£b can lead to further Rayleigh instabilities, depending 
on the parameters. Then, after reaching its maximum 
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Figure 20. Schematic view of the location of the differ- 
ent PE configuration types in the e/l^ 1 phase diagram. 
All simulations are performed at / = |. The dashed 
lines indicate the location of the simulation series B, C3 
and C4. 



extension at £b = 1.25a the chain slowly shrinks, since 
the enlarged Coulomb repulsion gets overcompensated 
by condensing counterfoils. Finally wc reach a collapsed 
conformation at £b = 9.0a. The non-monotonic behav- 
ior of the chain extension is qualitatively the same as 
in the good solvent case and qualitatively well under- 
stood ^01 CUD 02] ; however, the decrease is faster and 
more pronounced in the poor solvent case |19l I55j . 
Scaling theories have predicted that with the onset of 
condensation the pearl-necklace should collapse in a first- 
order transition |49ll^l lT0ll50| . However we find always a 
smooth distribution of countcrions which looks like a dis- 
tribution that can be calculated within PB theory |43| . 
The countcrions get pulled closer to the macroion as the 
Coulomb coupling is increased, and a rather high Bjer- 
rum length, or similarly, a large number of condensed 
counterions are needed to collapse the PE to a globule. 
Also the osmotic pressure does not show any dramatic 
decrease with £b, so we believe that for an adequate de- 
scription of the collapse we need a refined theory. 

A very interesting, and as we believe, new conforma- 
tional regime opens up basically for those £b values be- 
tween the maximal chain extension up to the collapsed 
state. Here we find conformations that are reminiscent 
of a cigar-like shape , but turn out to look more like a 
sausage for longer chains. Since with increasing £b more 
counterions are attracted towards the chain the pearl 
pearl repulsion is getting screened such that the pearls 
slowly coalesce and the conformation is not stretched on 
longer length scales anymore. Due to its shape we have 
termed it the sausage regime. Note that the actual con- 



formation depends also on / where smaller values of / 
lead to thicker sausages. Conformational snapshots are 
shown in figl21l 

At the crossing point of series B and C3 at / = |, 
£b = 1.5cr and e = 1.75&bT being in the pearl-necklace 
regime, we have studied the A m -dcpendence of the chain 
conformation. Starting from short chains which form a 
globule, e.g. one pearl, we enter a number of Rayleigh in- 
stabilities upon elongation. A few exemplary snapshots 
of series Al have already been shown in fig. 
As a next step we have studied the dependence of the 
conformations upon the charge fraction /. As one can 
see in fig. lathis opens up a completely new plane in the 
phase diagram of polyclcctrolytcs. The three simulation 
series CI, C2 and C3 all scan the entire ^B-range and 
differ only in the charge fraction /. 
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Figure 21. Schematic view of the location of the differ- 
ent PE configuration types in the £b/ f phase diagram 
The simulations are from series CI at / = 1, C2 at / = \ 
and C3 at / = \ performed at e = 1.75/cbT. 



In mean field theories for polyelectrolytes we find two 
important parameters. The first one is the Manning pa- 
rameter £ = tB^ x j which plays a role in all attempts 
to rcnormalize the charge of highly charged polyelec- 
trolytes 01] . The second one is a measure for the overall 
Coulombic repulsion on a charged chain €b& -1 / 2 - As a 
first Ansatz for a scaling theory for highly charged poly- 
electrolytes at finite density an effective charge fraction 
/ ff = j is often used, trying to combine both parame- 
ters. In this framework the three simulation series should 
behave identically. In contrast to this we find e.g. that 
the maximal extension i?B(max) for the three scries dif- 
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fers strongly. Namely we find for CI i?E(max) = 130<7, 
for C2 i?E(max) = 66.4cr and for C3 R^{max) = 13.5er. It 
is also striking that only the series with larger / has a 
regime where the chains are stretched and they behave 
as if they were in a good solvent. 

The Rayleigh instability occurs at the same value 
£b^ 1 ] 2 ~ 1/12 for all three series as expected because 
the counterions do not play a dominant role here. In- 
creasing £b from this point on the three series behave 
very differently. Whereas the series with f = h an d 
/ = 1 show a cascade of Rayleigh instabilities budding 
more pearls until they reach stretched conformations the 
series with f = \ has a dumbbell conformation at the 
maximum extension. The maximum extension itself is 
reached at different values for £b (sec fig. HOf) . For se- 
ries CI the maximum extension is probably more re- 
stricted by the chain entropy than that is determined 
by the interplay between repulsive and attractive inter- 
actions. The chains then slowly shrink and enter the 
sausage regime. Still at the same value of £b the chains 
with / = 1 are much longer than the chains with f = \ 
or / = -i. Finally the collapsed conformation is reached 
roughly at the same value of £b sa 8a. In this regime we 
have almost a dense electrolyte solution with mobile ions 
inside the chains so we should be close to the critical be- 
havior of a Coulomb fluid. This suggest that the collapse 
should occur at roughly the same value of the coupling 
parameter v m v c £s,/a which is the interaction energy of 
two oppositely charged ions at contact in units of k^T . 
In principle the phase space for polyelcctrolytes has far 
more than the shown three dimensions. As we know from 
previous studies ^3 EOI EH EH EH| also the density is a 
very relevant parameter and would have to be included 
into the phase diagram. Further important parameters 
are the valency of the counterions, and added salt con- 
centrations which we did not investigate at all here. This 
all reflects the fact that there is presumably no general 
parameter for the Coulomb interaction. 

8 Chain conformation and experiments 

The last section focuses on the connection between the 
chain conformation and experimentally accessible ob- 
servables like the characteristic ratios r and a and the 
form factor Si. 

8.1 Characteristic ratios 

The characteristic ratios r = (Re/Rq) 2 and a = Rq/R-r 
are often used as a first step to characterize conforma- 
tions that are extended in one dimension. Since both 
i?G and Rr are accessible with experimental methods it 
is interesting to know how a can be used to distinguish 
different conformation types. 

In fig. [221 an d fig- ESI we show the change of the charac- 
teristic ratios r and a with £b for the series CI, C2 and 
C3, e.g. for different /. For the globular conformation at 



C1:f=1.0 — 




2 4 6 8 10 



Figure 22. Change of the characteristic ratio r with 
£b for different / (scries CI, C2 and C3). 



£b = Oct we find r w 2.2 and a « 0.9. For dumbbells we 
find values for r between 4 and 4.5 and a « 1.6. With 
further increase of £b both observables reach a maxi- 
mum roughly at the maximal chain extension (compare 
to fig I10|> . Then r and a decrease monotonically until 
the collapsed conformation is reached at high values of 
£b where we find r ~ 3 and a ~ 1. For a globular object 
one can calculate r to be between 2 and 10/3 which is 
consistent with our findings. For a completely stretched 
object one would have r = 12 which is of course not 
reached with the simulations. The form of the shown 
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Figure 23. Change of the characteristic ratio a with 
Ib for different / (scries CI, C2 and C3). 

curves for r and a clearly indicates that it is not possi- 
ble to deduce the conformation type since at each value 
the chain could either be in a pearl-necklace conforma- 
tion or, at higher £b in a sausage like conformation. 
Since for the pearl-necklace regime itself a is a mono- 
tonic function of np we have plotted the values of a for 
all simulations in this paper that we assign to this regime 
in fig. 1241 One can see a clear jump of a between np = 1 
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Figure 24. Characteristic ratio a for pearl-necklace 
structures as function of np. The ellipses combine points 
close to np = 1,2 and 3. 



and np = 2. But the range of a between conformations 
with two and three pearls are already overlapping. With 
increasing np the slope of a(np) is decreasing and a fur- 
ther structure discrimination is not possible. So even in 
the pearl-necklace regime it is problematic to use a to 
discriminate pearl-necklace conformations with different 
numbers of pearls. 

8.2 Form factor 

More information about the chain conformation is con- 
tained in the form factor Si (see eq. |5J. In fig. l2~5l wc 
show Si(q) for different simulations from series C2. The 
figure shows the form factors for the different confor- 
mation types we have found in our simulations, namely 
the neutral globule, the dumbbell, the three pearl con- 
formation, stretched chains, sausage like conformations 
and a collapsed chain with most of the countcrions be- 
ing inside. The neutral globule shows Porod scattering. 
The strong oscillations show that the globule has sharp 
boundaries and does not fluctuate strongly. Si (q) for the 
dumbbell exhibits a shoulder at q ss O.ba^ 1 correspond- 
ing to a length of w 12.5<7 that can be identified as the 
distance between the two pearls. The strong decrease 
at higher q- values indicates again a Porod scattering in 
this case coming from the surface of the pearls. The 
minima of the Porod scattering are smeared out by fluc- 
tuations of the shape and size of the pearls. For the 
three pearl conformations we find two shoulders one at 
q w 0.2er -1 and one at q f» 0.4er _1 . They are also indi- 
cations of the pearl-pearl distances. But the shoulders 
here are already less pronounced than in the case of the 
dumbbell. The stretched conformation has a qualita- 
tive different Si(q) showing a scaling with g _1 over a 
large g-range. It follows that the chain is stretched on 
length scales larger than 10a up to its full length of 66a. 
For the sausage like conformation we can not identify 
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Figure 25. Form factors for the different conformation 
types while changing £b for series C2. The straight lines 
indicate a stretched chain with Si oc q^ 1 and Porod 
scattering with Si oc q~ 4 . 



any scaling regime or other signatures. Thus we can 
not obtain another length scale than the overall chain 
extension. The collapsed globule shows again, like the 
neutral globule, a strong decrease of the scattering in- 
tensity at large g-values. The rudimentary observable 
Porod scattering is strongly smeared out showing that 
the competition between attractive and repulsive forces 
induces large fluctuations on the surface of the globular 
object. 

After this overview over the scattering of the different 
conformation types we discuss the form factor of pearl- 
necklace conformations with a larger number of pearls 
in more detail. In this case the scattering is not only 
influenced by the fluctuations in shape and size of the 
pearls and strings but also by the fluctuations of the 
structure type, namely different numbers of pearls. 
The form factor shows four different regimes which re- 
assemble the different length scales that are involved in 
a pearl-necklace structure. For chains taken from series 
Al with N m = 382 the form factor is plotted in fig. EH In 
the Guinnier regime at Rcq •C 1 the radius of gyration 
Rq can be calculated from Si{q) =N m (l- (R G q) 2 /3). 
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Figure 26. Form factor Si for typical pearl- necklace 
conformations. The dotted lines are fits to certain q- 
ranges, see text. The marked region is enlarged in fig. EH 
The chains have a length N m = 382 and on average 
tip = 4.5 pearls. 



This yields Rq = 16.8a ± 0.3a in good agreement with 
the directly calculated value Rq = 16.9a ± OAa. In 
the following range 0.07(T _1 < q < 0.4<r _1 the single 
chain structure factor scales as Si(q) cx q~ x . The chain 
conformations are thus stretched on length scales larger 
than 15cr. At q f=a 0.5<t _1 one can see a weakly pro- 
nounced shoulder in Si. A closer look to this region 
reveals that S\{q) has an inflection point at q = 0.46(7~ 1 
which is shown in fig. 1271 A comparison with the in- 
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Figure 27. The pearl-pearl distance leaves only a very 
small signature in the form factor as can be seen in the 
close view of a comparison of the total form factor and 
the intra pearl form factor. 

tra pearl scattering shows that this is due to inter pearl 
scattering. This becomes more clear when one looks at 
the analytic scattering function S a of a linear arrange- 
ment of n homogenous spheres with a diameter rp and 
a distance rpp. S a is the product of the inter pearl scat- 



tering Sinter and the intra pearl scattering Si n t ra . The 
inter pearl scattering is given by 

e , oVV M sin(grppfc) 

Sinter = U + 2) [11 - k) : . (16) 

The intra pearl scattering is that of a homogenous sphere 
and thus given by the Porod scattering: 

_ sin(grp) - qr P cos^rp) 

"Jintra — 5 U'J 

grp 

Dividing out the intra pearl scattering from Si gives ac- 
cess to the inter pearl scattering and thus rpp. From 
the inflection point in the inter pearl scattering one can 
calculate the pearl pearl distance. This yields rpp = 
13.6cr, again in accord with the directly measured value 
?*pp = 13.3cr. In the high q-range between q ss 0.9(7 _1 
and q « 2.5cr~ 1 we find S\(q) cx q~ 4 , the typical Porod 
scattering with a small dip at q « 1.7tr _1 . Fitting the 
data to ea. H7l vields a pearl radius rp w 2.6a which again 
compares well to the directly calculated value rp ss 3.0cr. 
We conclude that the cooperative effect of fluctuations 
on overlapping length scales broadens all characteristic 
signatures which can be revealed by scattering under ex- 
perimental conditions like polydispersity and line charge 
density fluctuations. Thus necklaces might be difficult 
to detect. We find the most pronounced necklace signa- 
tures for the dumbbell conformation. 

9 Conclusion 

We have studied a dilute solution of strongly charged 
polyelectrolytes in a poor solvent by means of molecular 
dynamics simulations for a variety of different parame- 
ters. 

We have developed a cluster algorithm to characterize 
pearl-necklace structures even for single configurations 
and performed an extensive data analysis of all the pearl- 
necklace conformations in our simulations. 
We found that the range where scaling predictions are 
applicable is confined to a small range of parameters. 
One either has to stay close to the infinite dilution limit 
or use weakly charged chains to exclude countcrion ef- 
fects. The reason is that even in a dilute solution there 
is already a delicate interplay between the countcrion 
distribution and the chain conformation for strongly 
charged chains. This became in particular evident in 
our investigations studying the dependence on the sol- 
vent quality and the Coulomb parameters £b and /. 
Our results challenge approaches that apply Manning- 
like charge renormalization and simple scaling concepts 
to strongly charged polyelectrolytes. Since we found 
strong effects already for dilute solutions we expect an 
even more complicated countcrion chain coupling for 
semi-dilute and dense polyelectrolyte solutions. 
We have discussed the different types of fluctuations 
present for pearl-necklace structures and quantified them 
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for our data. We showed that the free energy differences 
between the structure types can be very small and that 
the size and position of the substructures exhibit large 
fluctuations. 

In a preliminary classification of the phase space of PEs 
in poor solvent we have shown that the phase space is 
multi-dimensional and that the region for pearl-necklace 
structures is rather small. In this context we have also 
investigated the Coulomb induced collapse of poor sol- 
vent PEs which appears to be a rather smooth transition. 
There we have shown that the strong screening in the 
vicinity of the charged chains leads to sausage like con- 
formations instead of pearl-necklaces. Our data suggest 
that the collapse is determined by the Bjerrum length 
since it occurs at £b ~ 8cr for all systems which lead us 
to suggest that in the collapsed state ion-ion correlations 
oc £b play the most prominent part. 
Finally we computed some experimentally accessible ob- 
servables like the form factor and characteristic size ra- 
tios. This should help to analyze the experimental data 
and to discover pearl-necklace signatures which, accord- 
ing to our understanding, are not very pronounced. 
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A Simulated Systems 

Here we give a detailed overview over the parameters 
of the simulated systems and some basic observables, 
namely the chain length N m , the end-to-end distance 
Re, the hydrodynamic radius i?u, the characteristic ra- 
tios r = (Rb/Rg) 2 an d a = Rg/Rh and the osmotic 
coefficient OC — Tl/pi g , where n is the pressure and pi g 
is the ideal gas pressure. If possible, also the number of 
pearls np and the average number of monomers in a pearl 
gp is given. In the tables to ED-Re, Rh and £b are given 
in a , e in kp,T . All other quantities in the tables are di- 
mensionless numbers. The statistical error for the chain 
extensions Re and i?H is smaller than 10%. gp exhibits 
a systematic error of ±4 monomers. For gp > 30 the 
number of pearls has an accuracy of ss 5%. For smaller 
gp values this error becomes larger. The statistical error 
of the pressure calculation gives rise to an error of the 
osmotic coefficient of ±0.06 resulting in a large relative 
error for small values of OC. 



TABLE 1: Series Al: N m = 48. ..478, 

e = 1.75fc B T, £b = l-5cr, / = 1/3, p c = 1.0 x lO^V" 3 



N m 


Re 


7? H 


r 


a 


np 


gp 


OC 


48 


3.48 


2.19 


3.25 


0.88 


1.00 


46.0 


0.96 


94 


4.56 


2.71 


3.40 


0.91 


1.01 


93.9 


0.68 


142 


12.9 


4.11 


4.63 


1.46 


1.97 


68.8 


0.60 


190 


15.0 


4.60 


4.49 


1.54 


2.01 


90.4 


0.44 


238 


25.7 


5.86 


5.84 


1.83 


2.94 


75.5 


0.46 


286 


30.4 


6.51 


6.10 


1.89 


3.24 


82.5 


0.42 


334 


38.2 


7.44 


6.73 


1.98 


3.94 


78.6 


0.40 


382 


45.4 


8.24 


7.24 


2.05 


4.53 


78.0 


0.35 


430 


55.4 


9.22 


7.55 


2.19 


5.22 


75.2 


0.36 


478 


59.7 


9.81 


7.61 


2.20 


5.60 


78.3 


0.39 
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TABLE 2: Series A2: iV m = 100.. .300, 

e = 1.75&BT, In = 1.5cr, / = 1/2, pc = 6.7 x lO^a" 3 



N m 


Re 


i? H 


r 


a 


np 


gp 


OC 


99 


26.8 


5.59 


7.17 


1.79 


2.77 


23.3 


0.54 


199 


68.1 


10.43 


8.20 


2.28 


4.13 


20.7 


0.32 


299 


102.9 


14.37 


8.88 


2.40 


5.99 


18.9 


0.39 



TABLE 3: Series A3: iV m = 100.. .300, 

e = 1.75fe B T, 4s = l-5cr, / = 1/2, p c = 6.7 x lO" 5 ^ 3 



iVm 


R E 


R H 


r 


a 


np 


gp 


OC 


99 


22.1 


5.08 


6.72 


1.68 


2.89 


25.1 


0.46 


199 


54.9 


9.15 


8.24 


2.09 


5.25 


21.4 


0.36 


299 


83.4 


12.6 


8.69 


2.24 


7.68 


20.4 


0.35 



TABLE 6: Series C2: £ B = 0.0...6.0er, 

/ = 1/2, e = 1.75/cbT, N m = 199, p c = 5.0 x lO^cr" 3 



in 


R E 


R B 


r 


a 


np 


gp 


OC 


0.0 


4.41 


3.27 


2.22 


0.91 


1.0 


199 


0.97 


0.167 


4.63 


3.30 


2.39 


0.91 


1.0 


199 


0.93 


0.333 


14.2 


4.59 


3.72 


1.60 


1.87 


102 


0.82 


0.5 


26.6 


5.72 


5.38 


2.00 


2.90 


61.4 


0.79 


1.0 


54.7 


8.90 


7.53 


2.24 


4.80 


25.0 


0.57 


1.5 


65.2 


10.4 


8.22 


2.19 


— 


— 


0.45 


2.0 


66.4 


10.8 


8.40 


2.12 


— 


— 


0.36 


2.5 


62.6 


10.5 


8.52 


2.04 






0.27 


3.0 


59.9 


10.3 


8.48 


2.00 






0.21 


3.5 


54.2 


9.84 


8.26 


1.92 






0.16 


4.0 


48.0 


9.23 


8.08 


1.83 






0.13 


6.0 


27.8 


6.96 


7.03 


1.50 






0.06 


9.0 


7.45 


4.37 


3.16 


0.96 


1.14 


173 


0.05 



TABLE 4: Series B: e = 0.0k B T...2.0k B T, 



N m = 238, £ B 


= 1.5(7 


/ = 


1/3 , Pc 


= 5.0 x 10" 5 C7- 3 


e 


Re 


R B 


r 


a 


np 


gp 


OC 


0.00 


102 


16.2 


8.88 


2.11 






0.66 


0.50 


100 


14.9 


9.28 


2.21 






0.62 


1.00 


82.2 


13.1 


8.81 


2.11 






0.63 


1.25 


68.7 


11.0 


8.44 


2.15 






0.53 


1.35 


58.5 


9.78 


7.91 


2.13 


5.70 


25.2 


0.52 


1.45 


45.9 


8.22 


7.32 


2.07 


5.12 


35.6 


0.53 


1.55 


34.3 


6.89 


6.56 


1.94 


3.92 


52.9 


0.45 


1.65 


24.2 


5.83 


5.76 


1.73 


2.91 


76.8 


0.48 


1.75 


16.8 


5.03 


4.72 


1.54 


2.14 


108 


0.36 


1.85 


13.8 


4.68 


4.36 


1.41 


1.90 


123 


0.31 


2.00 


11.2 


4.34 


4.21 


1.26 


1.52 


155 


0.26 



TABLE 7: Series C3: £ B = 0.0.. .9.0(7, 

/ = 1/3, e = 1.75fc B T, N m = 199, p c = 5.0 x lO^V- 3 



£b 


Re 


i? H 


r 


a 


np 


gp 


OC 


0.0 


4.41 


3.27 


2.22 


0.91 


1.0 


199 


0.97 


0.25 


4.59 


3.29 


2.37 


0.91 


1.0 


199 


0.93 


0.5 


5.10 


3.31 


2.86 


0.91 


1.0 


199 


0.76 


0.75 


10.5 


4.04 


4.01 


1.30 


1.62 


120 


0.64 


1.0 


13.0 


4.39 


4.24 


1.44 


1.87 


104 


0.54 


1.25 


13.5 


4.51 


4.22 


1.46 


1.98 


97.4 


0.47 


1.5 


13.5 


4.50 


4.38 


1.43 


1.95 


99.2 


0.38 


2.25 


12.2 


4.47 


4.32 


1.31 


1.57 


125 


0.29 


3.0 


12.3 


4.44 


4.90 


1.25 


1.30 


152 


0.22 


4.5 


9.86 


4.29 


3.99 


1.25 


1.06 


188 


0.10 


6.0 


7.77 


4.06 


3.39 


1.04 


1.01 


198 


0.15 


9.0 


6.35 


3.82 


3.07 


0.95 


1.00 


199 


0.07 



TABLE 8: 



Series C4: £ B = 1.5. ..12.0(7, 



TABLE 5: Series CI: 

/ = 1, e = 1.75fc B T, N m = 200, 



£ B = 0.0.. .12.0a, 

Pc 



£b 


Re 


Rh 


r 


Q 


np 


gp 


OC 


0.0 


4.41 


3.27 


2.22 


0.91 


1.0 


200 


0.97 


0.083 


15.8 


4.70 


4.04 


1.67 


2.0 


95.3 


0.94 


0.25 


72.1 


10.6 


7.99 


2.41 


3.83 


21.1 


0.86 


0.5 


111.3 


14.9 


9.51 


2.42 






0.75 


0.75 


124.8 


16.3 


9.88 


2.44 






0.63 


1.0 


129.5 


17.0 


9.98 


2.41 






0.55 


1.25 


125.8 


17.1 


9.70 


2.36 






0.46 


1.5 


130.2 


17.1 


10.2 


2.39 






0.40 


3.0 


104.8 


15.9 


9.34 


2.16 






0.18 


6.0 


67.66 


13.0 


8.22 


1.82 






0.06 


8.0 


31.94 


8.46 


6.74 


1.45 






0.06 


10.0 


9.204 


5.45 


3.04 


0.97 


1.0 


200 


0.04 



= 1/3, e = 


l.5k B T, N m 


= 94, 


Pc = 


1.0 x 


lO-V- 


£b 


Re 


Rh 


r 


a 


np 


gp 


OC 


1.5 


13.5 


3.93 


5.33 


1.49 


2.04 


41.2 


0.68 


1.8 


13.9 


4.00 


5.47 


1.48 


2.10 


39.6 


0.62 


2.1 


15.1 


4.15 


5.88 


1.50 


2.24 


36.2 


0.59 


2.4 


15.7 


4.23 


6.14 


1.50 


2.27 


35.3 


0.49 


2.7 


15.6 


4.22 


6.41 


1.46 


2.26 


35.9 


0.53 


3.0 


15.6 


4.23 


6.47 


1.45 


2.26 


35.8 


0.40 


3.3 


14.4 


4.10 


6.43 


1.39 


2.17 


38.5 


0.45 


3.6 


13.2 


3.97 


6.22 


1.33 


2.06 


41.6 


0.30 


3.9 


12.2 


3.85 


6.17 


1.28 


1.92 


45.6 


0.33 


5.1 


8.6 


3.44 


5.08 


1.10 


1.38 


66.7 


0.22 


5.4 


8.2 


3.39 


5.06 


1.08 


1.30 


71.2 


0.17 


5.7 


8.0 


3.37 


4.96 


1.07 


1.26 


73.6 


0.08 


6.0 


7.7 


3.34 


4.81 


1.05 


1.22 


76.4 


0.12 


9.0 


5.7 


3.14 


3.58 


0.97 


1.03 


91.4 


0.19 


12.0 


5.1 


3.06 


3.17 


0.94 


1.01 


93.4 


0.04 
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TABLE 9: Series D: £ B b - 1 f 2 = 0.25 =const., 
e = l.75k B T, N m = 200, p c = 5.0 x IQ-^a^ 



iB f 


Re 


i? H 


r 


a 


np 




OC 


0.25 


1.00 


59.4 


9.54 


7.81 


2.23 


4.88 


21.6 


0.75 


1.00 


0.50 


32.3 


6.66 


6.96 


1.84 


4.37 


39.2 


0.44 


2.25 


0.33 


8.3 


4.15 


2.90 


1.18 


1.11 


179 


0.17 


4.00 


0.25 


5.9 


3.55 


3.00 


0.96 


1.00 


201 


0.13 
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